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In this paper we analyze the influence of transport mechanisms (diffusion and sedimentation) 
on the structure of monolayers of particles irreversibly adsorbed on a line. We focus our attention 
on the dependence of the radial distribution function g(r) and the saturation coverage Oao on the 
gravitational Peclet number N g . First, we study the probability density of adsorption onto an 
available interval using approximate solutions of the transport equation and computer simulations. 
Combining our results with an approximate general formalism, we can obtain values of #00 and the 
^~ ) ", gap density, which agree with our simulations. We also show that, for large gravity, the coverage 

■ #00 approaches the ballistic limit following a power law in N g that is independent of the number of 

I dimensions, as has been observed in simulations pl[ . 
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' I. INTRODUCTION 

The adsorption of particles of colloidal size (macromolecules, latexes, bacteria, etc.) from fluid suspensions to solid 
surfaces is a complex phenomenon of great interest. Much effort has been devoted to the study of the mechanisms 
involved in this process pj. A complete description should consider the different steps involved: i) transport of the 
particles from the bulk suspension toward the interface; ii) interaction with the substrate (including the layer formed 
by the previously adsorbed particles); iii) surface diffusion and desorption. For many colloidal particles, for example 
some proteins |2||| or latexes neither surface diffusion nor desorption is observed in the time scale accessible 

to experiments: the particles remain immobilized after adsorption, and the process can be considered irreversible. 
Consequently, non-equilibrium configurations are generated and, when the surface coverage attains a given value, a 
1 jamming configuration is obtained, with no space available on the surface for the adsorption of new particles. 

The simplest model that attempts to describe these irreversible adsorption phenomena, is the Random Sequential 
1 ' Adsorption (RS A) model |^,[^,^|j9| Jl^Jll| , p^| . In this model, particles are sequentially added to the surface by iteration 
of the following algorithm: a) a random position is selected for the addition of a new particle; b) if the new particle 
overlaps any particle already adsorbed on the surface, the adsorption attempt is rejected; c) if the new particle does 
not overlap any other particle, then the adsorption attempt is accepted; d) the time necessary for the adsorption of the 
particle is proportional to the number of adsorption attempts. The RSA model has been extensively studied in both 
discrete and continuous surfaces ||. Exact results for the kinetics, jamming coverage and distribution functions can 
be obtained analytically for one-dimensional surfaces |Tc|| ; for two-dimensional surfaces one has to use approximate 
j_j ■ methods pd| , p!2[ or computer simulations ||. 

Clearly, the RSA model does not consider the transport of the particles toward the surface. Instead it assumes 
that new particles arrive to the neighborhood of the surface with uniform probability, and then they interact with 
the partially covered surface on the basis of excluded volume interactions. More realistic models must consider the 
different transport mechanisms that are present in colloidal systems: diffusion, gravity, externally imposed flows, and 
hydrodynamic and double layer forces JL3[ . 

For large enough particles suspended in a fluid at rest, gravity is the dominant force: large particles sediment 
following a ballistic trajectory toward the surface. This situation is described by the ballistic deposition (BD) model 
|i~4|.[l5]], in which step b) of the RSA model is modified in the following way: if the incoming particle overlaps a 
pre-adsorbed one, then it rolls down the steepest descent path until either it reaches the surface and is adsorbed, or it 
is trapped in a cavity over other particles, and is rejected. This model is exactly solvable for one-dimensional surfaces 
SLpJl , and extensive approximate studies have been realized for two-dimensional surfaces [14[ . 
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In the opposite limit of very small particles, Brownian motion becomes the dominant transport mechanism. A model 
describing this situation is the diffusion RSA model (DRS A) JItJ : the initial position of each particle is randomly 
chosen in a plane at a given distance of the surface, and its trajectory is simulated by using a Brownian dynamic 
algorithm |l8|; if the particle touches the surface it adsorbs irreversibly. Furthermore, to avoid unbound Brownian 
trajectories, particles arriving at points too far from the surface are rejected. The simulation of this model is much 
more expensive computationally than the RSA model and, furthermore, it admits no exact solution even for one- 
dimensional surfaces. Surprisingly, the jammed state obtained is very similar to the jammed state obtained with the 
RSA model: very small differences appear in one-dimensional surfaces, while in two-dimensional surfaces the jamming 
states obtained from both models are indistinguishable with the precision of the simulations. Analytically, good 
approximations can be obtained for jammed one-dimensional surfaces from approximate solutions of the diffusion 
equation |19|. 

For particles of general size both, the Brownian motion and gravity have to be considered. Simulation studies have 
been done including the gravity force in the Brownian dynamics algorithm pCfl . The results show how the properties 
characterizing the surface structure change smoothly with increasing gravity, from the DRSA behavior in the limit of 
small gravity to the ballistic values in the limit of large gravity. One interesting result is that the dependence of the 
jamming coverage on the size of the particles seems to scale to a common curve both for one and two-dimensional 
surfaces |f2l|j . Therefore, the information obtained from one-dimensional surfaces can be relevant also to the more 
realistic two-dimensional ones. 

Finally, any realistic model has to include the effect of the hydrodynamic interactions (HI) , which become important 
when the particles are in the neighborhood of the surface. Simulation results [^2 23 suggest that the effect of HI on 
global averaged quantities is small although they should be taken into account for a fine analysis of the local structure. 
It can be useful therefore to study simpler models that neglect hydrodynamic interactions but can be analytically 
solved. 

Our aim in this paper is to develop an analytic approximation to the description of the adsorption of particles on one- 
dimensional surfaces considering the gravity force and Brownian motion and neglecting hydrodynamic interactions. 
To do this, we first need solutions for the transport equation in the presence of a flat adsorbing surface and many 
particles fixed on it. These solutions can be obtained with a sufficient approximation when only one particle is 
fixed on the surface. Then, a superposition approximation has to be done to study surfaces at finite coverage; this 
superposition has shown good results for the DRSA model Jl9[ , and in the present model it also gives good agreement 
with the simulations. 

In section || we describe the model and the analytical tools we use in detail. In section III the one-particle effects 
are studied by approximate solution of the transport equation and computer simulation. The results of this section 
allow us to obtain an approximate description of the adsorbed phase that is compared with simulations in section[F\~ 
Finally, the appendices A and B show how the approximate solutions of the transport equation can be obtained. 



II. DESCRIPTION OF THE MODEL 



We want to study a simple model that may allow us to understand the effect on the structure of the adsorbed layer 
of two simple transport mechanisms, namely diffusion and sedimentation. We consider an adsorbing surface at z = 
and a semi-infinite fluid in the region z > 0. Spherical particles of radius R suspended in the fluid diffuse in the XZ 
plane and sediment due to the effect of a uniform gravitational field in the Z direction. If the center of a new particle 
arrives at the z = line, it is adsorbed irreversibly at the contact point. We assume that the bulk concentration of 
particles is so small that interactions between them are negligible. Consequently, each particle adsorbs independently 
of the other particles in the suspension, and the process can be considered as sequential. However the concentration 
becomes large at the surface, where adsorbed particles accumulate and interact via excluded volume effects with 
incoming particles. 

This adsorption problem can be studied in two steps: first, the transport problem for a single particle near the 
surface in the presence of the previously adsorbed particles has to be solved; then, the evolution of the adsorbed phase 
can be studied from the obtained adsorption probabilities. 

Let P(r,t) be the probability density to find the center-of-mass of a particle at the point r at time t; under the 
conditions of negligible inertia and small relaxation time, the transport of the particle to the adsorbing line is governed 
by the Smoluchowski equation, 



dP(r,t) _ ~ 



DVP(r,t) + uP(r,t)z , (1) 



dt 

u being the sedimentation velocity and D the diffusion coefficient assumed constant. The probability flux is given by 
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J(f, t) = -L>VP(r, t) - uP(r, t)z. 



(2) 



The hard-particle interaction between bulk particles and adsorbed ones is considered in the boundary conditions 
for eq. (|l|) by assuming that the radial probability flux must vanish at the exclusion surface (figure [l]) delimited 
by the preadsorbed particles. At the adsorbing line we impose perfectly adsorbing boundary conditions (P = 0). 
Furthermore, we assume the initial condition P(f, t = 0) = S(z — zq), zq being the initial distance of the center of the 
particle to the line; the value of z is not important provided that Zq > 2R. 

Equation (Q) with its boundary conditions can be made dimensionless by measuring distances in units of the 
diameter 2R of the particles and time in units of the characteristic diffusion time Tdiff — 4R 2 /D. Thus, the solution 
of (Q) depends on a single dimensionless parameter, the gravitational Peclet number N g , defined as the quotient 
between Tdiff and the characteristic sedimentation time Tdet = 2R/u, 



_ I d l£1 _2Ru_ SirR+gAp 
9 ~ r det D 3k B T ' [6) 



where g is the acceleration of gravity, Ap the difference between the densities of the particles and the fluid, and T the 
absolute temperature. If N g ^> 1 the motion of the particles is deterministic and we recover the Ballistic Deposition 
model (denoted in the following by BD), whereas if N g <C 1 Brownian motion predominates and we recover the DRSA 
model |1S[ |. Note that N g is proportional to R 4 , and therefore we can define a dimensionless radius as 



the factor 2 has been introduced to agree with the definition given in [J20 21 1. For polystyrene beads in water at 
T = 300K, Ap = A5>kgjrr? and therefore R* ~ 0.817i? if R is expressed in microns. We will use these dimensionless 
quantities in the following. 

It is not possible to analytically solve the diffusion problem ([!]) in the general situation, when several adsorbed 
particles are present. We will limit ourselves to obtentaining approximate solutions in the presence of only one 
adsorbed particle. To describe the general situation, then, we will introduce a superposition approximation by 
assuming that the adsorption probability at a given point depends only on the position of the nearest particles, that 
is, those limiting the free interval of the line to which the point belongs. One expects this hypothesis to fail only for 
small gaps (and therefore only slightly affecting the layer structure) if gravity is not too high. If gravity is strong as 
compared to diffusion, when an incoming particle touches an adsorbed one, it is expected to adsorb in its immediate 
vicinity and the superposition hypothesis holds also for the smallest gaps. 

According to the previous discussion, the solution of the transport equation will lead to an adsorption model in 
which particles are sequentially deposited onto an infinite line with an adsorption probability density that depends 
on the gap in which adsorption takes place. This means that the adsorbing particle only significantly interacts with 
the particles that are limiting this gap, and interactions with other particles are not considered. The solution of 
the Smoluchowski equation (fil) gives the adsorption probability at any point of an available interval. The problem 
of characterizing the structure of the adsorbed monolayer can be studied starting from the kinetic equation for the 
evolution of the free gaps. We summarize the method below, and a more detailed discussion can be found in ]l9| . 

Let G(h, t) be the number density of gaps with length h at time t and k(hf, h) the probability per unit length and 
per unit time that the deposition of a particle in a gap of length h! > 1 produces gaps of length h and h! — h — 1. The 
governing kinetic equation for the irreversible adsorption process is [ fl9| (see figure [l]) : 

=-k (h)G{h,t) + 2 f dh'G(h' t t)k(h',h), (5) 

"t Jh+1 

where ko(h) is the total rate at which gaps of length h are destroyed by the addition of a new particle: 

k (h) = / dtik{h 7 ti). (6) 
Jo 

The gap distribution in the jamming state can be obtained without knowing the detailed time evolution. From the 
balance equation (|^) it is possible to derive a time-independent equation for the total number density of gaps of length 
h that have been created at any time along the process, n(h) [|19| : 

poo 

n{h) = 2 dtip(h',h)n{ti), (7) 

Jh+l 
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where p(h', h) is the probability density that the first particle arriving at an interval of length h' creates two new free 
intervals of length h and h' — h — 1, 



From the definition of k (h), eq. (^J), this function satisfies the normalization 

rh'-l 



p(h',h)dh = l. (9) 
In addition eq. (0) must be supplemented by a normalization condition for n(h). From the definition of n(h) we have 

n 

1 

(1 + h)n(h)dh = 1. (10) 

A global measure of the particle packing is the coverage defined as the relative length (area in 2 dimensions) of a 
line of total length L covered by N particles of radius R: 



2RN 



(11) 



The coverage increases monotonically (due to the irreversible nature of the process) until a saturation value 9^ , when 
no available space remains for adsorption of new particles. In this situation only gaps with h < 1 survive and the 
saturation coverage can be obtained from n(h) using the fact that the particle number is equal to the number of gaps, 

000= f n(h)dh. (12) 
Jo 

Let us remark that eq. (fjj) expresses the fact that the total number of gaps of length h can be computed from the 
number of gaps with length h! > h + 1 and the probability p(h', h) of obtaining an interval of length h from an interval 
of length h! . In principle, from eqs. (|7|)-(|l2"|) it is possible to obtain n(h) and the saturation coverage 9oo if p{h' , h) is 
known. For each specific adsorption model p(h' , h) must be previously known. 

In order to obtain p(h', h), we must solve eq. (ED) with two adsorbed particles limiting a gap of length h! according 
to the hypothesis assumed in the derivation of eqs. The accuracy of this ansatz is analyzed by means of 



Brownian dynamics simulations in section IV 



III. ADSORPTION PROBABILITY 



In this section, we study the probability density of the adsorption of a new particle on a line in which one adsorbed 
sphere is already present. To find this probability we compute approximate solutions of the Smoluchowski equation 
(Q) and we perform Brownian dynamics simulations. The results obtained will be used in the next section to determine 
the structure and coverage of jammed lines. 

To use the integral equation (f?j), we need the probability density p(h',h) defined in (@), which can in principle be 
obtained by solving the Smoluchowski equation (|l|) in presence of two preadsorbed particles separated by a center 
to center distance h' + 1. However, the solution to this problem is very difficult and therefore we have adopted 
a superposition hypothesis assuming that the joint effect of the two preadsorbed particles can be obtained from a 
superposition of the one-particle effects. 

The one-particle adsorption problem is defined as follows: the center of an adsorbed particle is fixed at the point 
(x, z) — (0, 0) and, at time t = 0, a new particle starts to move from an initial point at a given distance from the 
surface zq, with a random initial value of x. The initial height zq is immaterial and it will be assumed to be infinite. 
In the presence of a non-vanishing gravity field, the particle will be adsorbed with probability one. We normalize 
P(r, t) in such a way that the initial probability per unit length is equal to one. Therefore, the probability density of 
adsorption at a point of the adsorbing line z = is: 

/>oo 

p(x)=- J z (x,z = 0,t)dt, (13) 

J 
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where J z (x, z, t) is the z component of the probability flux defined in eq. (|J). Using the superposition hypothesis, we 
write for the probability density of adsorption on a gap of length h! : 

p(h',h) ~N(h')p{h + l)p(ti -h), (14) 

where N(h') is determined by the normalization condition (^J). To obtain p(x) we do not need to solve the time- 
dependent Smoluchowski equation ([!]). We only need the time-integrated probability density, 

/>oo 

i/)(f)= / dtP(f,t), (15) 



V^(r) + N g ^ = 0. (16) 



that obeys the equation (obtained by the integration of (Q) and using dimensionless units) 

dz 

The boundary conditions satisfied by ip are the perfect adsorbing boundary condition on the line z=0: 

i/)(x, z = 0) = 0, (17) 
and the condition that the total probability flux coming from z = oo is 1 per unit length, 

d ^ 1 N g rj}{r) -> 1, z -» oo. (18) 



dz 

Using eq. ( ^3|) and the boundary condition (|l7|), p(x) can be obtained from ijj as 

c?-0 



az 



(19) 



z=0. 



Another boundary condition is necessary to account for the presence of the previously adsorbed particle in the origin 
of coordinates. The hard sphere interaction between diffusing and adsorbed particles implies no radial flux at the 
exclusion sphere (of radius 2R) centered at the origin; therefore we have (see figure |l|): 



dil){r,9) 



+ N g ij){r = l,0)sin0 = O. (20) 

r=l 



dr 

The solution ip(r) can be split into two parts, 

i/;(f) =^(°)(f) + ^ (1) (r), (21) 

where ip^ is the solution corresponding to a clean adsorption line, and tp^ reflects the presence of the adsorbed 
particle and must vanish at large distances. The solution of ([11]) for ip^{f) with conditions ( p^Jl^ ) is 

^)( XlZ ) = ±( 1 - e - N s*). (22) 
The probability density of adsorption in the presence of the fixed particle is obtained using eq. ([l9j): 



p(x) 



dz 



+ 



dz 



= l + pW, (23) 

z=0 



where p^ is the deviation of p from the uniform distribution. In the interior of the region excluded by the fixed 
particle (|x| < 1), the total flux has to be zero and therefore p^ = —1. In the available region (|x| > 1), p^ gives the 
local increase of the adsorption probability due to the particles that have been in contact with the adsorbed particle; 
these particles are rejected in the RSA model, but here are allowed to diffuse and will adsorb in the neighborhood of 
the fixed particle. The presence of the fixed particle does not change the total adsorption probability, so the integral 
of p^ over the entire line has to vanish. Therefore one has 



P 



Wdx = 1. (24) 



The equation and boundary conditions for ip^(r) can be obtained by substitution of (|2l|),(p2|) in (|i6[) , jT7|) , (p0[) . It 
is not possible to obtain a solution for general values of N g , and different approaches have to be used in the limits of 
small and large gravity. 
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A. Small gravity 



We obtain a multipolar expansion for p{x) by solving the equation for ip^(r) and substitution in (|23|). To obtain 
ipW(r) we use the Green function method, for the detailed calculations see Appendix A. The (formally exact) result 
for p(x) can be expressed in the form 



p(x) = 1 + — - 
ir\x\ 



oo 

E 

n=l 



(2n - l)tp2n-l 



A' 



\x\ > 1, 



(25) 



where K n (x) is the Bessel function of third kind and the coefficients <02n—i satisfy the infinite linear system of equations 

■El3B- _ . _ 

By substitution of the multipolar expression (E5h in jl4[ ) we obtain the probability density of adsorption on a gap 
of length h! for given N g . The normalization factor N(h') is obtained by numerical integration of (^). For practical 
uses, the series (E5j) must be truncated by using a finite number of terms; the larger the gravitational Peclet number 
N g is, the larger the number of terms that have to be retained in (p5f) to achieve a given precision. We recall that 



if we want to obtain the saturation coverage Ooc we have to solve the integral equation (H) with the kernel p(h', h) 
given by eqs (|l4]) , (|25|) . If N g is of order 1 we need only a few terms in the summation and the numerical solution of 
the integral equation is possible, but for intermediate to high values of N g we need a very large number of terms and 
this expansion becomes useless. Therefore, it is necessary to obtain an asymptotic expression for large values of N g 
to enable a description in the entire range of variation of N g . We will return to this question in the next subsection. 

We can see how the DRSA result [[L9) can be recovered by taking the N g — > limit in the previous expression. For 
small N g we need only one term in (|25|), because f>2n-i — § <^in + 0{N g ). Therefore, 



p{x) 



+ 0{N g ) wl + l/x 2 



(26) 



The last approximation has been obtained using the small argument expansion of the Bessel function K\{x) ~ -. It 
gives the known zero-gravity result |]l9| and we see now that it gives an approximation for small gravity valid in the 
region i < 2/jV 9 . 

The first moment of p^ gives the mean distance to the origin where particles touching the fixed one adsorb. From 
eq. (p5|), and using the well-known formula for the derivatives of Bessel functions, K n -x{y) + K n+ \{y) = —2K' n (y) 
we obtain: 



< x >i = 



(27) 



For N g —> only the first term of this sum has to be retained, obtaining 



2 K (N g /2) 

Ng Kl{Ng/2) 



log 



N B 



(28) 



The mean distance increases without bound, as can be expected from the DRSA limit given by eq. (pq). 



B. Large gravity 

For large values of N g one can try to solve eq. (|l^) by a perturbative approach. In the N g —> oo limit the higher 
order derivatives in eq. ( |l6|) are negligible, and an approximate equation is obtained retaining only the advective term. 
This is equivalent to completely neglecting diffusion. The corresponding solution satisfying the boundary condition 
at infinity is ip = 1/N g , which does not verify the boundary conditions at z — and r — 1. The reason is that, no 
matter how large N g is, diffusion becomes important near the boundaries of the system: ip must change appreciably 
in a thin boundary layer, and the corresponding derivatives are no longer negligible in that layer. 

The problem can be solved using singular perturbation techniques as discussed in appendix B. The result is that, 
for large values of N g , p^\x) scales to a function of the form 

p^{x)^Ny^[{x-l)Nl% (29) 



6 



where function cj) 1 - 1 ' is defined in the appendix. The mean distance at which a particle that hits the preadsorbed one 

2 /3 

can be adsorbed is therefore a quantity of order N g , 

f°° A 
< \x\ >i = / p {1 \x)xdx w 1 + 



/•OO 



/•OO 

/ <£ (1) 1-106.... (30) 

J o 



Obviously this distance decreases when iV s grows because diffusion is weak and the particles adsorb closer. If N g = oo, 
a particle hitting the fixed sphere will roll over it and be adsorbed in its immediate vicinity (at |x| = 1) as corresponds 
to the rolling mechanism of BD. If N g is large but finite, diffusion disturbs the deterministic motion and the new 

2 /3 

particle can be adsorbed at some distance from the adsorbed one, leavingbetween them a gap of size h ~ N g ' . 

We can give a simple reasoning supporting the scaling given by eq. ( p0|) for large values of N g . Particles falling 
over the fixed particle will roll down over it and diffuse away along a thin layer of thickness 6. This thickness can be 
estimated from the condition that the radial probability flux due to gravity, which is proportional to N g sin has to 
be compensated with the diffusion flux that can be estimated as 5 . Therefore, in dimensionless units the boundary 
layer thickness is 

N g sm 

This thickness increases when the particle approaches the line 9 = 0. For small angles the surface of contact between 
the particles becomes vertical, and the rolling picture is no longer applicable. Instead, we can imagine that particles 
diffuse away from the fixed sphere at a certain angle, 9q. This separation is expected to occur when the x coordinate 
of the center of the falling particle is beyond the region excluded by the fixed one. Given eq. (|3l|) , this implies 

[1 + S(9q)] cos 9q — 1. Solving this equation for small angles, we obtain 0q ~ (2/Ng) 1 ^ 3 . From this point the incident 

1/3 

particle falls a distance zo — (2/-/V s ) ' ; the dimensionless time needed to reach the z = line is to ~ z o/N g = 

2 1 / 3 N g 4 ^ 3 . Due to the horizontal Brownian motion performed during the particles' falling, one expects a horizontal 
deviation of order 5x ~ \/2to — (2/N g ) 2 ^ 3 . This simple argument gives an exponent and order of magnitude that 
agrees with the more exact singular perturbation result (Ejof). 



C. Brownian Dynamics Simulations 



Computer simulation allows us to study the adsorption probability for the whole range of N g values and to verify 
the validity of the results obtained in the previous subsections. The algorithm developed here will be also used in the 
following section to study the structure of jammed lines. 

We consider, as described in section a suspension of Brownian particles that is dilute enough to assume sequential 
adsorption. One by one, the center-of-mass of a new particles starts its motion at a fixed vertical position zq and 
at a horizontal position xq chosen at random from a uniform distribution. The value of zq does not affect the final 
adsorption probability as far as in the initial position there is no possibility of interaction with the adsorbed particles, 
that is, zo > 1. In our simulations we take zq = 1 to minimize the simulation time required. 

The particle motion is discretized in time in the following way. At every time step At the particle develops two 
motions performed sequentially in the simulations. It travels a vertical distance Azdet = —N g At (in the dimensionless 
units introduced in section |J) as corresponds to the sedimentation under gravity and a stochastic displacement Ar rc j 
as corresponds to the Brownian motion. Therefore, 

Ar = -N g Atz + Ar rd . (32) 

As is well known, the Brownian displacement in two dimensions follows a normal distribution law with zero mean and 
variance 

< Ar r 2 d >= 4Ai. (33) 

After these two motions are generated, the simulation time is increased by At and the algorithm is repeated until the 
particle reaches the adsorbing line z = 0. The time step At is chosen in such a way that ensures a typical displacement 
of the particle to be of magnitude e which is small compared to the particle diameter, e 1. The typical Brownian 
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displacement is of order V At, so At < e . On the other hand, the deterministic displacement is —N g At, so At < e/N g . 
Therefore, we impose 

At = min [e 2 ,e/N g ] . (34) 

The value of e is restricted by computer time limitations. All the simulations reported in this paper were obtained 
with e = 0.05. 

To complete the simulation algorithm, we need a collision rule to use in the case that an incoming particle touches 
the pre-adsorbed one. We consider two collision rules, depending on whether at the collision instant the particle is 
performing a deterministic or a random displacement. If the collision occurs during the deterministic displacement, 
we consider that the incoming particle rolls over the preadsorbed one. The collision rule for the Brownian motion is 
a perfect reflection of the particle trajectory, as in an elastic collision with an infinitely massive particle. This rule is 
an approximation, valid if the Brownian displacement is sufficiently small in relation to the radius of the boundary. 

When the center of the particle reaches the z = line, we consider that the particle adsorbs at the contact point, 
which is recorded. By iteration of this algorithm one obtains a histogram of adsorption frequencies as a function of 
the distance to the fixed particle. 

For small values of N g the results obtained by Brownian dynamics simulations agree with the analytical results 
obtained from the superposition approximation eq. (|l4|) , and the multipolar expansion cq. (|2^) , see figure |2|. For large 
values of N g the simulated distribution of adsorbed particles approaches the scaling (^) predicted by the asymptotic 
approximation, see figure ^. 

To obtain a simple analytic approximation in the intermediate to high N g regime, we fit the simulation results for 
p(x) with an exponential distribution, 

p(x) ~ l + be-^-V, \x\>l. (35) 

Although the asymptotic form of the adsorption probability obtained in appendix B is not exactly an exponential, 
this is a simple fitting function. To determine the value of b(N g ) from the simulation results we impose the condition 
that the thickness of the exponential has to be equal to the thickness of /jW obtained from the simulations. This is 
the condition of maximum likelihood for an exponential distribution [p8| , and it implies 

- =< \x\ >i -1. (36) 

For large values of N g we expect, based on our asymptotic argument (|30|), a dependence of the form < \x\ >i= 
1 + a N g .To determine the constant a we computed the mean distance of the adsorbed particles to the previously 
adsorbed one for different N g numbers. The simulations were performed for N g = 60, 80, 100, 140, 200, 400, 600 using 
40, 000 particles in each simulation and we obtained: 



N, 



2/3 
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b=^—, a = 1.08..., (37) 
a 

in agreement with the value 1.106... obtained from the asymptotic solution (|30|). In figure [| we compare this 
exponential fit, eqs. ([35j), (|37]) with the results of the Brownian dynamic simulations for N g = 60. 

Making use of the superposition hypothesis ( |l4| ) and the normalization condition (JoJ) the probability density of 
adsorption in a gap can be expressed in the simple form: 

p(h', h) = —^1 + ° e — - . (38) 

^ i h i + 1 _ 2e -b(h'-l) ^ > 

Although the exponential distribution does not reproduce all the details of the simulated p[x), it approaches its 
main characteristics, and the corresponding kernel ( ^8| ) is simple enough to allow a numeric evaluation of the integral 
equation (||). 

IV. STRUCTURE OF THE ADSORBED LAYER 

In this section we present results concerning the structure of the ads orbed layer at jamming, characterized by the 



coverage #oo and the radial distribution function g(r). In subsection IV A we study the variation of the jamming 
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coverage with N g and in subsection IVB we analyze the local structure of jammed surfaces described by g(r). The 
results obtained by numerical integration of eq. @) using ( j38|) or are compared with simulation results. 

The characteristics of the simulations are as follows. For each value of N g we have filled 1000 lines of length 
L = 100 and we have computed the mean number of particles adsorbed at saturation. The statistical uncertainty on 
the coverage (error bars) is estimated by a 95% confidence interval. In principle, particle trajectories are simulated 
as long as the surfaces are not entirely covered, but there exist some special cases that demand more attention to 
avoid prohibitive computer time. Consider a particle diffusing in the region between two pre-adsorbed ones, being 
the interval of the line between these pre-adsorbed particles of length less than 1. This particle only can adsorb if 
it obtains sufficient thermal energy to overcome the gravitational force and the geometrical barrier delimited by the 
pre-adsorbed particles. The time needed to escape from these "geometrical barriers" grows rapidly with N g . For 
N g ~ 1 the particle can return to the bulk after a small time and try to adsorb elsewhere. However, for sufficiently 
high N g this time becomes very large (infinite in the case of BD). In the practical situation, these particles play no 
role because other particles, coming from the bulk, will adsorb in the free gaps. There are different ways in which 
this effect can be taken into account |2q] . In our simulations we adopt the following criteria: we choose a maximum 
residence time r and, if a given particle remains time t > r in a trap, it is eliminated from our simulations and a new 
particle starts the motion following the previous rules. We take a different characteristic residence time r for each 
N g . For small N g only a few particles are eliminated, but for N g ^> 1 usually all particles that fall into intervals of 
length h < 1 are eliminated. For R* < 1 we take r = 15, r = 8 for 1 < R* < 2 and r = 1 for R* > 2. We have 
checked that our results are not sensible to increments of r. Finally, the last stage is simulated in the following way. 
When only gaps with length hi < 3 survive, the initial position x n is chosen at random between the centers of the 
particles limiting these gaps, at z = 1 as for the other particles. When only targets of length h! < 2 remain (at most 
one particle can be adsorbed in each) we fill these by using a RSA algorithm, i.e. by randomly placing a particle in 
it without any transport process. Only a few particles are deposited following these two last rules and therefore the 
structure (characterized by g(r)) is not significantly altered. 



A. Coverage at the jamming limit 



For adsorption on an infinite line, the saturation coverage can only be a function of the unique dimensionless 
parameter present in the transport equation, namely N g . It is more convenient to represent the variation of 9^ as a 
function of the dimensionless radius, R*, defined in eq. (|4|), which gives a more compact scale than N g . 

In figure ^| we show the variation of 9^ with R* obtained by simulation (crosses) and by solving the integral 
equation (^) following the numerical method described in p9[ . The solid line in the region R* < 1.5 (N g < 10) has 
been obtained with the kernel given by substitution of the multipolar expression ( p5|) truncated to five terms in eq. 
(|l4|). We also solve, for all values of R*, the integral equation (Q) with the kernel obtained using the exponential 
approximation, eq. (|3|); the result is represented by the solid line for R* > 1.5 and by the short dashed line for 
R* < 1.5. 

Three main regions can be recognized in the graph: 

1. R* < 1: In this region the Brownian motion dominates. The coverage varies slowly with R* and its value 
remains close the DRSA value (9^ ~ 0.7506). For example, for the R* = 1 case we obtain 9oo ~ 0.7539 by 
means of simulations and the numerical solution of the integral equation with the multipolar kernel yields 
6oo ~ 0.7528. Simulation results agree with the curve obtained using the multipolar solution of transport 
equation. The curve obtained using the exponential kernel ( |38| ) gives smaller values for 9^ and in the limit 
of vanishing gravity it approaches the RSA coverage 9^ ~ 0.74759 . . . instead of the DRSA value. 

2. 1 < R* < 2.5: 9oo rises quickly with R*. In this region the numerical solutions of the integral equation 
with the multipolar expression and with the exponential approximation ( |38| ) are close. However, they show 
a slight discrepancy with the simulation results that are systematically greater. 

3. R* > 2.5 6 X approaches slowly to the BD limit (9^ D = 0.80865...), for example, if we take R* = 10 
(N g = 10 4 ) simulations give 6*00 = 0.807.... We obtain a good agreement between the simulations and the 
solution of (0) using the exponential kernel (^q) . A simple expression describes the asymptotic approach 
to the strong gravity limit, N g — > 00: 

MA 9 »1)~ 0.2445 Ng- 2 '\ (39) 

This expression is shown in fig. [| by dashed lines, and can be obtained considering a simplified model. 
For large values of N g , the explicit form of seems not to be of great importance provided that it 
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gives a value for the mean distance between the adsorbed particles and the pre-adsorbed one in agreement 
with eq.(|3^). We can use instead of the exact density, an approximation in which is described by a 
normalized step function of thickness 2 < |a;| >i= 2AN g 2 ^ 3 in agreement with 
solvable 0|, giving (|3^) in the limit N g — > oo. 

We have performed some control simulations to study the origin of the slight discrepancy between the solid line and 
simulations in region 2 on the graph. This effect is due to the fact that for this range of N g the probability density 
of adsorption on a gap not only depends on the gap's limiting particles but can also be influenced by third neighbors. 
This can be clearly shown by performing simulations with three pre-adsorbed spheres and comparing the resulting 
adsorption density probability with the result of the superposition approximation (|lj). In figure g we show the 
adsorption probability density in a situation with three pre-adsorbed spheres on a line of length L = 20 and R* = 1. 
Between spheres one and two there is a gap of length hi — 2.5 and between spheres two and three there is a small 
gap of hi = 0.05 that does not allow particle adsorption therein. Simulation results show an asymmetric adsorption 
probability, which clearly displays the effect of the third particle. The solution corresponding to the superposition 
hypothesis ( |l4|) , which assumes that the adsorption probability in one gap is given by the product of the one particle 
solutions corresponding to the particles at the ends of the interval, shown in the continuous line, is symmetric. We 
also note that the simulation results fit well with a distribution constructed as a product of the three one-particle 
solutions (dashed line). For larger R* values, the ansatz of independence of gaps holds, and the theoretical expression 
shows no significant deviation from their simulation counterpart. 



30 ) . This model is exactly 



B. Radial distribution function at the jamming limit 

The radial distribution function, g(r), reflects the correlation existing between the adsorbed particles. The definition 
adopted here is the usual one in Statistical Mechanics (see for instance |26|]). In this subsection, we present histograms 
of g(r) in the jamming limit, obtained from the Brownian dynamic simulations used in the previous subsection to 
determine #oo- 

In fig.0 and |s| we show g(r) for various values of R* . At the jamming limit, g(r) presents a logarithmic singularity 
at r — 1, like the RSA model. The simulation results, however, are histograms that have a finite value at contact 
representing an average of g(r) over a finite distance interval. This contact value is nearly constant for R* < 1 whereas 
it grows with R* for R* > 1. 

For R* < 1 our g(r) are indistinguishable from the results obtained for DRSA (R* = 0). The g{r) for R* = 1 is 
shown in fig. |t], note the disappearance of correlations between particles after a few diameters. We recall that the 
variation of with R* in this region is less than 1%, as we pointed out in the previous section. Therefore, we can 
conclude that for R* < 1 the structure of the adsorbed layer is almost independent of gravity. 

In fig. | we show the radial distribution function for N g = 30 (R* = 1.968) and for N g = 240 (R* = 3.310). A 
comparison between figures and |^ shows that peaks in g(r) increase and are steeper as R* increases, revealing the 
tendency of large particles to pack closer together than smaller ones. This is because diffusion is a small effect for 
large particles, and there is an increased probability of adsorption close to already adsorbed particles. We can observe 
that the smaller the effect of gravity (as measured by R*), the poorer the structure in the ra dial d istribution function. 



For 1 < R* < 2.5 not only the coverage changes rapidly with R* (as noted in subsection IV A) but g(r) does also. 
Note for example that, for R* = 1.968, two secondary peaks are visible whereas only one is visible for R* = 1. For 
R* > 2.5 the obtained g(r) are very close to the correlation function of the ballistic deposition model, although the 
delta function singularities present in the BD model at r = 1, 2, . . . are smoothed due to the effect of diffusion. Our 
model differs from other generalized models [jl6| in which the delta funciton singularities appear only for R* = oo. 

At the jamming limit g(r) cannot be obtained completely using the formalism presented in section but it can 
be obtained for 1 < r < 2 noting that n(h) is equal to g(h + 1) for < h < 1 except by the normalization. By the 
definition of these functions one has 

n{h)=e 2 00 g{h+l), h<l. (40) 

The numeric solution of eq.(Q) performed in the previous subsection to find 9oo also gives the gap density n(h). The 
result agrees with the radial distribution function obtained by simulation according to eq. , as shown in fig. ^ for 
R* = 1 and in fig. §for R* = 1.968, 3.310. 
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V. CONCLUSIONS 



A 1+1 dimensional adsorption model has been analyzed to investigate the influence of transport mechanisms 
(diffusion and sedimentation) on the structure of irreversibly adsorbed monolayers of particles. We have studied a 
simplified model in which hard spheres suspended in a two-dimensional fluid adsorb onto a line. The results may be 
relevant to understand the more realistic case of adsorption on a surface, if the suggested scaling behavior between 
the 1+1 and 2+1 dimensional cases holds |^lf . An argument in favor of this hypothesis for large gravity is presented 
below. 

By means of Brownian dynamic simulations we have obtained the saturation coverage and the radial distribution 
function g(r) for different values of the gravity number N g . These results have been compared with those obtained 
from an approximate analytic formalism: first we analyze the probability density of adsorption onto a gap with the 
help of the transport equation and Brownian dynamic simulations; once we know this probability distribution, the 
saturation coverage 9oo(R*) and the density of gaps at the jamming limit can be obtained numerically by using an 
integral equation developed in a previous paper |l9[ . 

We have introduced a superposition hypothesis, assuming that the adsorption probability at a given point depends 
only on the limiting particles of the gap in which the particle adsorbs. As showed by simulations, this hypothesis fails 
when the gap is small and 1 < R* < 2, but this affects doo only slightly. 

The probability density of adsorption in the presence of one particle was obtained by solving the transport equation 
as a multipolar expansion. This expansion is useful for low to intermediate N g values, but not for high N g numbers. In 
the large N g regime we have shown, using singular perturbation techniques, that the probability density of adsorption 
p^\x) has the scaling form (|29|). The simulations show that this scaling is approximately satisfied even for moderate 
values of N g . Instead of the exact scaling function, we have useda simple exponential function for the numerical 
computations that fits the simulation results. The corresponding numeric solution of the integral equation (Q) agree 
well with the simulations even for R* > 1. 

The properties of the adsorbed layers at jamming are as follows. For R* < 1, oo (R*) varies less than 1% and g(r) 
is close to the DRSA (R* = 0) form. Therefore, the effect of gravity is nearly negligible in this regime. For R* > 1, 
oo (R*) grows quickly with N g and g(r) displays peaks that emerge more significantly. The peaks in g(r) increase and 
are steeper when R* grows reflecting the tendency of large particles to pack closer than smaller ones. For R* > 2.5 
the saturation coverage slowly approaches the BD value following the asymptotic expression (|3^). In this regime g(r) 
is very close to its BD counterpart, although in our model g(r) takes finite values at r — 1, 2, . . . due to the effect of 
diffusion instead of the delta function singularities of the BD case. 

It can be easily verified that the same perturbation techniques used in appendix B can be applied to the 2 + 1 
dimensional case, leading to the scaling given in eq. (p9l) with the same function This fact provides a physical 

basis, at least for large gravity, for the scaling of 9oo\R*) in a common curve for the 2 + 1 and 1 + 1 cases. In the 
asymptotic limit the relevant property of the adsorption probability is its thickness, and one can replace p^ 1 ' with a 

2/3 

step function with the same thickness, A = kN g . We expect that the jamming coverage approaches the ballistic 
limit linearly in A, and therefore 

«=^a-i-<«^, (4i) 

where ay is a constant that depends on the dimension of the system. Thus, with an adequate change of scale the 
curves of oo (N g )/9^ D in 1+1 and 2+1 dimensions are coincident, at least for large values of N g . This scaling has 
been observed in simulations and seems to apply for all values of N g [^lj . 

In a more realistic extension of our model, hydrodynamic interactions must be considered in the transport equation, 
but we expect that their effect is small. In the DRSA case, hydrodynamic interactions make the adsorption 
probability nearly uniform due to the enhanced mobility parallel to the surface and therefore the coverage is expected 
to be close to the RSA case. In the case of BD it has been shown that while the jamming coverage does not 
change significantly, the local structure is strongly affected by the hydrodynamic interactions. 
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APPENDIX A: MULTIPOLAR SOLUTION OF TRANSPORT EQUATION 

In this appendix we obtain ip^{r) using the Green function method. The equation and boundary conditions for 
■0W(f) can be obtained by substitution of in @,@,(||). If we write ^(r) = /(r)e _JV »*/ 2 , the equation 

for f(r) is 



with boundary conditions: 



Or 



N 2 

v 2 /« = -fm, 



f{x,z = Q) = a 

+ ^sm8f(r = 1,9) = -sinf 



The Green function corresponding to (Al) and satisfying the boundary condition ( ]A2[ ), G(r, r'), is defined by 

iV 2 

V 2 G(r,f') - -JLG(r,f") = 5(f-f') 
G{x,z = 0,r') = 0. 



(Al) 

(A2) 
(A3) 



(A4) 
(A5) 



Using the method of images, it is easy to obtain the solution to this equation, that can be expressed using Bessel 
functions of the third kind: 



K a [-2-\?-?'\ \-K [ -^-\f-r 



N a 



(A6) 



r = (x 1 z) r ={x,z) r* = (x,—z). 
We can obtain an integral equation for f(r) making use of the properties of the Green function. If we multiply eq. 



(Al) by G(f,f') and eq. (A4) by —f(r) and add the results, we obtain: 

G(f, r ')V 2 /(r) - /(r)V 2 G(r, r') - -f{f)5{r- r'). 



(A7) 



Now, integrating (A7) over the whole volume V, limited by the adsorbing line at z — 0, the excluded volume of the 
pre-adsorbed particle and a closing surface S, and then using the Stokes theorem, we obtain: 



f(?) = - dV [G(r,r')V 2 /(f) - /(r)V 2 G(f, r ')] 



dS 



G(r,r')Vf(f) 



dS 



/(r)VG(r,f) 



(A8) 



If the surface S goes to infinity and we use the boundary conditions for / and G(r, r'), the only contribution to the 
surface integrals comes from the excluded area of the preadsorbed particle. Therefore, using dS = Rd9 = d6 we 
obtain: 



fir') 



dO 



r' > 1. 



(A9) 



r=l 



To take profit of this expression, we use the Fourier expansion of the Green function (A£) that can be obtained 
using a well-known theorem for the Bessel functions |E9| : 



K o ( Pyii + r 2 - 2rir 2 cos0 J = K m (Pr 2 )I m (/3ri) cosm^, r 2 > 



n ■ 



m— — oo 



Use of (A10) in yields: 



r I sm(n9) sin(n#'), r < r 



(A10) 



(AH) 
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Substitution of this expression in eq. (A9) yields 



f(f')= — ^K n (^r')I n (^)sm(ne') / d0sin(«0)/ 



TV ^ 2 ' A 



N N 

I n+l{ — ) + In-l(-) 



sin(n0') / d0f{r = l,0)sin(r 
Jo 



(A12) 



Now, we note that our problem is symmetric under the change x — > tt — x, and, therefore, only the integrals with odd 
n contribute. In add itio n, if we apply the boundary condition of null radial flux at the excluded area defined by the 
preadsorbed sphere (A2) and define ip(8) = f(r = 1,0) we obtain 



/(*") 



- y^K 2n -i(— r')sin(2n- 1)0' / <20sin(2n-l) 
n ~, 2 Jo 



n=l 
N 



sin We 2 



U ( N 3) + 



If we define the Fourier-sine transforms: 

92n-l(N g ) 



dO sin(2n- 1)6 sin 6e— su 

r 2n _i(iV ff ) = / d6sm(2n-l)6smdip(6), 
Jo 

V2n-x{N a ) = [ d9sm{2n - 1)0^(0), 



finally one obtains for f(r, 6) (we omit ' for simplicity) 

2 ~ 

/(r,0) = - > ip 2n -i 

7T » — ' 



2^, K 2n .^r) , 

^ — sin 2n — 1)0. 

AWi(^) 



(A13) 

(A14) 
(A15) 
(A16) 

(A17) 



n=l 



By substitution of r = 1 in (A13) and making use of definitions (A14-A16) we obtain the following equations for the 
coefficients <f2n—i' 



V2n-1 



T2n-1 



9ln-\ 



N„ 



z z 1 - -it K ^-\\^r)\h n -2\^r) + hni-f)) 



-4(2m-l)(2n- 1) 



7T (2n + 2m - 3)(2n + 2m - l)(2n - 2m + l)(2n - 2m - 1) 



(A18) 
(A19) 



The probability density for the adsorption of a particle can be obtained in a form of a multipolar expansion by using 
eq.flAl7l) with the definition of f{r), V> (1) (V) = /(f) e^* 2 / 2 in definitions @, @: 



P\x) = 1 + -j-, > (2n - l)<p 2 »_i— 7^— 



b > 1. 



(A20) 



For small values of N g , one has <72n-i = f^ni + C(JV ff ), and all the terms in eq.(A20) except the first vanish to order 
0(JV fl ), giving ^ 2 „_! = f <5 nl + C(A 9 ). 



APPENDIX B: BOUNDARY LAYER SOLUTION OF THE TRANSPORT EQUATION 

Our aim in this appendix is to obtain an approximate solution of the transport equation ( |lt| ) valid for large values 
of N g . When N g is very large, the laplacian term in the transport equation becomes negligible in the entire domain 
except near the boundaries, where the value of the derivatives can be large. We have then a singular perturbation 
problem ]3C| ] . An approximate solution in this limit can be reached by matching different approximations in the bulk 
and in the boundary layers that can appear near the boundaries. 



13 



The first approximation to the outer solution can be obtained neglecting the laplacian term in the transport 
equation. The resulting equation is simply dip/dz = 0, and the solution satisfying the boundary condition at infinity 
eq. ([l?]) is a constant, 

VW = l/N g . (Bl) 

Note that this is an exact solution of the transport equation, and therefore the expansion of the outer solution ends 
here. 

This outer solution does not match the boundary conditions at z = and r = 1, and therefore some boundary 
layers have to be introduced. 

(I) . A first boundary layer, appears near the line z = 0, where the value of ip has to change from the outer value, 
l/Ng, to the value imposed by the boundary condition, 0. Derivatives with respect to z become large, and the term 
d 2 ip/dz 2 in the laplacian is no longer negligible. The boundary layer thickness, e, can be estimated from the condition 
that the derivatives with respect to z are the dominant terms in the transport equation and have the same order. 
This implies e = 0(N~ r ). Therefore, the first approximation to the inner solution in this region is the solution of 

d 2 ^ 1 dtb 1 

^ +Ar £Jm =0 . (B2) 
dz z dz 

The solution satisfying the boundary condition at z — and matching the outer solution for z 3> N^ 1 is 

^n = ^-(l-e-^ 2 ). (B3) 

This is also an exact solution of the transport equation, so the expansion of the inner solution in this boundary layer 
ends here. In the absence of any adsorbed particle, no additional boundary layers are present and, in fact, we have 
reproduced the exact solution for an empty line given by eq. (^2|). 

(II) . In the presence of one particle in r = 0, a new boundary layer appears near the bo und ary of that particle, 



r = 1. The function rj) has to change in a small boundary of thickness 8 from the outer form (Bl), in which there is a 
radial probability flux, to a form in which that flux vanishes at r = 1 . The order of magnitude of the different terms 
of the transport equation can be easily evaluated by writing it in polar coordinates, and considering that each radial 
derivative introduces a factor 0(S~ 1 ): 

Ify 1 d 2 V . cV costfcty 

H 7T- + -5-^7 o + N q sm9— + N q — =0. 

dr r dr r 2 80 2 9 dr a r dO 

0{5- 2 ) 0(5- 1 ) 0(1) OiNgS- 1 ) 0{N g ) (B4) 

The first and fourth terms of this equation are the largest when N g — > oo, and the appropriate distinguished limit is 
obtained when they have the same order of magnitude, that is, when 6 = 0{N g ' 1 ). To obtain the solution in this 
boundary layer, we define the scaled inner variable n = N g (r — 1), and we expand the function ip*£ in successive 
approximations, ifj 1 ^ = ijj^ + + • • •. The first approximation to the inner solution is then, retaining the terms of 
order 0(N 2 ), 

sin 9-j—=0. B5 

drf drj 

The boundary conditions for cause the normal flux to vanish at r\ = and the solution has to match the outer 
solution for r] —> oo. The first integral of eq. (pq) is immediate, 



^—+sm6^ =C(9). (B6) 
dn 

The term on the right-hand side is an integration constant that, in principle, can be an arbitrary function of 9. 
However, the left hand side is proportional to the radial flux and, by the boundary condition at r] — 0, the integration 
constant must vanish. We have, then, a homogeneous first-order equation for ip^\ whose general solution is 

V> (1) = A(<9)e-" sine . (B7) 

It is not possible to obtain the "integration constant" A{0) by direct matching with the outer solution, because this 
first approximation to the inner solution vanishes when rj — > oo. As we will see, the reason for this is that the first 
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approximation t/>W is a quantity of order 0(1), whereas the outer solution is ©(TV" 1 ). Matching is possible only for 
the second-order approximation, i/^ 2 \ which is 0{N~ 1 ). Retaining terms of order 0{N g ) in the transport equation 
we obtain the equation satisfied by 

1 suiff— — = —N„ — hcosf 



drj 2 drj 9 \ drj 89 

= -N^cosOe- 11 ^ 9 \A'{9) - A(9)(r)cos9 + tan0)] . (B8) 

The first integral satisfying the boundary condition at rj = is 

<9i/> (2 



drj 



+ sin^ (2) = -Ng 1 cose / dr]' e - ri sin 9 [A' (d) - A{d){r]' cos 6 + tan 9)} . (B9) 

Jo 



It is not necessary to go beyond this calculation to obtain the matching condition with the outer solution. It suffices 
to recognize that the left-hand side of the last equation is —Ng 1 times the radial probability flux, and this quantity 



has to match when 77 — s- cxd with the outer value for this flux, which from (Bl) is simply — sin#. We obtain the 
condition: 

n (A'(0)-A(6)tzn.6 A(9)cos9} 

-sia9 = cos9<—^ — L (BIO 

\ sin(9 sin 2 9 J V ' 

This is a differential equation for the unknown function A(9), whose general solution is A{9) = sin# + £?tan#. The 
arbitrary constant B has to vanish in order to obtain a finite probability density for 9 = ir/2. This completely 
determines the first approximation to tp(^. 

However, this approximate solution is not valid for all values of the angle 9. The reason is that the boundary layer 



thickness is S = (N g sm9) , as is evident from eq. (B7), and it increases with decreasing angle. Therefore, the 



boundary layer approximation bre aks down when 9^0. It is easy to verify that the order of magnitude of the last 



term on the left-hand side of eq.(B4), which has been neglected, becomes comparable to the retained terms when 
9 = 0(N g ±/J ). Therefore, the range of validity of the solution obtained is 

= sin0e-" sine + OiN- 1 ), 9 » N^ 1 ' 3 . (Bll) 

-1/3 



(III). For values of the angle 9 = 0{N g 1 ), the last term in eq. (|B4|) has to be included in the boundary layer 



equation together with the two terms retained in the previous approximation. All these terms have the same order 
of magnitude when the radial thickness of the boundary layer is S = 0(N g 2 ^ 3 ). These range of values for the angle 
and the radius define a new boundary layer; we introduce new appropriately scaled variables, £ = (r — l)iV 2 ^ 3 , and 
t = 9Ng^ 3 . We introduce also an asymptotic expansion of the probability density, ipf^ 1 — fii(N g )(f)( 1 ' > +/i2(iV g )^^ 2 - ) + • • •, 
where the coefficients fj,i(N g ) are infinitesimals of increasing order. The equation for <fP-> is, retaining the terms of 
larger order in N g : 



(i) 



t-j— = 0. (B12) 



dr d£ 2 d£ 

We now have a parabolic equation that has the form of a time-reversed diffusion equation with a drift proportional 
to the "time" r. For r — > 00, the solution has to match with the previous boundary layer solution ipf£, valid for large 
values of the angle. It is easy to see that this implies 

H 1 (N g ) = N- 1 / 3 , 

-> re~ r? , r -► +00. (B13) 
That is, for large values of r we recover a normalized "Boltzmann distribution" , corresponding to neglecting the 



r-derivative in eq. (B12). This condition, together with the boundary of condition of zero normal flux at £ = 0, 
d(j>^/dt; + T(f>^ = 0, completely specifies the function (j)^. Furthermore, the condition of zero flux at £ = 0, implies 



that the integral of <fP-> with respect to £ is a conserved quantity, and as a consequence of eq. (BIS) it has value 1: 



0«(£,rK = l. (B14) 



(i 
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We have not been able to analytically solve eq. ( B12| ). However, it is easy to obtain a numeric solution satisfying 
the cited boundary conditions. The solution is a positive normalized function, vanishing exponentially for £ — > oo, 
and with a well-defined limit when r — > 0, c^ 1 -* (£, r = 0) (see figure |]). An important quantity is its first moment, 

A=/ O °°dee^(0 = 1-106.... _ 

The adsorbing boundary condition at z — 0, however, is still not satisfied by this approximation. From the previous 



— 1/3 

(£, 0) and this result, obtained from an approximation valid for 9 — 0{N g ), 



solution we have lim T ^o — N, 

has to be modified to agree with the boundary condition ip(8 = 0) = 0. Indeed, when 9 — 0(N g ' 1 ), we are in the same 
situation as in the boundary layer (/): the derivatives with respect to 9 are of order e = 0(N^ 1 ), and the boundary 
layer equation is eq. (B2). The range of validity of the boundary layer approximation ip 1 ^ 1 is therefore 

tin =N- 1/3( t> {1) (£,T), N- 1 <9<1. (B15) 



(IV). A new boundary layer is needed for 9 = 0(N 1 ), in which the approximate transport equation is eq. jB2j ) 



The solution of that equation satisfying the boundary condition at 6 = is 



IV 



C(0 (1 - e- N °°) . 



(B16) 



The arbitrary function C (£) is obtained by imposing matching with the approximate solution in the previous boundary 
layer ipin^ ■ We obtain 



^(r,^)^iV- 1 /3^)(e,r = 0)(l 



-N a 6 



(B17) 



This function gives the correction to the probability flux, originated from the presence of the fixed particle, p 



(i). 



P {1) {r) 



dz 



W(e,r = 0), £ = (r-l)iV 2 / 3 . 



(B18) 



z=0 



Therefore, the numerical solution of eq. (B12) specifies the final probability distribution. One important conclusion 

is th at, fo r large values of N g , scales to a function of the form Ng<f>^ [(r — l)Ng^ 3 ]. It is easy to verify, using 
eq. (B14), that this distribution is normalized, 



p (1 \r)dr 



0(i)(£ j7 - = O)d£ = l. 



(B19) 
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FIG. 1. Illustration of the governing equation for the non-uniform deposition 

FIG. 2. Probability density of adsorption in presence of a pre-adsorbed particle as a function of distance for N g = 10. Solid 
line: multipolar solution of transport equation (25) including 5 terms. Dots: Brownian dynamic simulations. 

FIG. 3. Probability density of adsorption in presence of a pre-adsorbed particle as a function of distance for N g — 60. Solid 
line: exponential distribution (35). Dots: Brownian dynamic simulations. 

FIG. 4. Logarithmic of the scaled probability density of adsorption in the presence of a pre-adsorbed particle. Simulation 
results for N g — 60 (crosses) and N g — 300 (triangles). Continuous line: Numerical solution of boundary layer equation. 

FIG. 5. Comparison between the saturation coverage Ooo as a function of R* obtained with several methods, (a) Dots: 
Brownian dynamic simulation, (b) Solid line: numerical solution of the integral equation (7) using the multipolar solution (25) 
truncated to five for R* < 1.5 and the exponential fitting (35) for R* > 1.5. (c) Dashed line: asymptotic expression (39). (d) 
Short dashed line: continuation of the numeric solution with exponential fitting (35) for R* < 1.5. 

FIG. 6. Probability density of adsorption in a line of length 40i? in the presence of three adsorbed particles, at positions 
xi — 13R, X2 = 20R, X3 — 22.1R. (a) Dots: Brownian dynamic simulations, (b) Solid line: superposition of the two 
one-particle solutions corresponding to the particles limiting each gap. (c)Dashed line: Superposition of the three one-particle 
solutions. 

FIG. 7. Radial distribution function g(r) for R* = 1. Dots: Brownian dynamic simulations. Solid line: Calculation of g(r) 
for 1 < r < 2 using the numerical solution of the integral equation (7) for n(h). 

FIG. 8. Simulation results for g(r) corresponding to R* — 1.968 (N g = 30), squares and R* = 3.310 (N g = 240) displaced 4 
units, crosses. Solid line: Calculation of g(r) for 1 < r < 2 and N g = 30 using the numerical solution of the integral equation 
(7). Dashed line: the same for N g = 240. 
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